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^ ■ ABSTRACT 

Context. The formation of intracluster light and of the extended halos around brightest cluster galaxies is closely related to morpho- 
3 ■ logical transformation, tidal stripping, and disruption of galaxies in clusters. 

Aims. Here we look for observational evidence to characterize these processes, by studying the morphology and kinematics of the 
£ — ' diffuse light in the core of the Hydra I cluster. 

Methods. We analyze Ks- and V-band surface photometry as well as deep long-slit spectra, and establish a link between the structures 

in the light distribution, the absorption line kinematics, and the LOS velocity distributions of nearby galaxies and planetary nebulae 
i 1' (PNs). 

Results. The central galaxy NGC 331 1 is surrounded by an extended symmetric outer halo with n=10 and an additional, off-centered 
envelope at R ~ 50" to the North-East. Its luminosity L v = 1.2 x 10'° (±6.0 x 10 s ) L Q corresponds to ~ 50% of the luminosity of the 
symmetric halo in the same region (~ 15% of its entire luminosity). The velocity dispersion of the halo rises to cluster core values, 
~ 400 - 500 km s~ , for R > 20". Based on measured PN velocities, at least part of the off-centered envelope consists of high-velocity 
accreted stars. 

We have also discovered two tidal streams in the cluster center, emerging from the dwarf galaxy HCC 026 and from the SO galaxy 
HCC 007. The HCC 026 stream is redshifted by ~ 1200 km s _I with respect to NGC 3311 (V N 33ii - 3800 km s" 1 ), similarly as 
HCC 026 itself, a fraction of PNs in the off-centered envelope, and several other dwarf galaxies nearby. The stars in one of the 
HCC 026 tails are known to be consistent with the low-metallicity population of HCC 026, and our photometry shows that this galaxy 
is already mostly dissolved in the tidal field. 

The tidal stream around HCC 007 extends over at least ~ 110 kpc. It is fairly thick and is brighter on the side of the asymmetric 
.< ■ outer halo of NGC 331 1, which it may join. Its luminosity is several 10 9 L Q , similar to the luminosity of the stripped-down galaxy 

HCC 007. The redshift of the stream is determined from a few PN velocities and is similar to that for HCC 007 and HCC 026. 
JZ. \ Conclusions. An entire group of small galaxies is currently falling through the core of the Hydra I cluster and have already been 

partially dissolved by the strong tidal field. Their light is being added to the outer halo and intracluster light around the cD galaxy 

NGC 331 1. The Hydra I cluster provides a vivid example of morphological transformation and tidal dissolution of galaxies in clusters. 

Key words. galaxies:clusters:general - galaxies:clusters:individual (Hydra I) - galaxies:kinematics and dynamics - galax- 
^ ies:individual (NGC 3311) 

o ' 

(N ; 

,— 1 , 1 . Introduction such as tidal interaction and harassm ent, mergers and cannibal- 

ism, gas stripping and starvation; see Poggianti ( 20041) : lDe Lucial 



X 



Galaxy clusters are the most massive virialized structures in the , 2 007) forreviews. Which of these mechanisms takes the leading 
universe, and may consist of thousands of galaxies. One of the role for a given galaxy morphological type in different environ- 



s' . most interesting fields in modern cosmology is the study of the mental cond itions still remains to be understood. 
?3 mechanisms for the growth and evolution of such systems and 

the evolution of galaxies within them. The hierarchical model In the nearby universe, the evolution of clusters as a whole 
predicts that structure formation and evolution occurs through and the evolution of galaxies in clusters can be addressed by 
the merging of smaller units into larger systems, and this model stu dying the dynamics of the intracluster light (ICL). The ICL is 
has been supported by many observational evidences. The obser- the diffuse light in galaxy clusters emitted by stars which are 
vational appearance of clusters and their galaxies additionally not boun d to any specific galaxies; for a review on the sub- 
depends on the evolution of the baryonic component, which is j ec t see lArnaboldi & Gerhard! d^OlOb . Wide field surface pho- 
less well understood. As galaxies fall into dense environments, tome try shows structures in the IC L on all scales, from a few 
their evolution is affected by a variety of dynamical processes arcminutes t o deg r ees on the sky (iThuan & Kormendvl fl977t 

iMihos et alj 120051; iRudick et all. l2009h . Recent studies have 

Send offprint requests to: M. Arnaboldi, e-mail: marnabol@eso.org shown that the ICL provides direct evidence for the dynamical 

* Based on observations collected at the European Organization for status of galaxy cluster c ores (IGerhard et al.l2007HDohertv et all 

Astronomical Research in the Southern Hemisphere, Chile, under the 120091; IVentimiglia et alj 1201 ll) . because it contains the fossil 

observing programs 082.A-0255(A), 076.B-0641(A), 065.N-0459 record of past interactions, due to its long dynamical time. 
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Cosmological hydro-dynamical simulations predict that the 
ICL is formed by stars that are unbound from galaxies during 
the interactions they experience as they fall through the cluster 
potential well and interact with other cluster galaxies and the 
cluster tidal field. In these simulations the ICL shows significant 
subst ructures on all scales in its spatial and velocity distribu- 
tion (iNapolitano et al.ll2003tlMurante et al.li2004tlWillman et alJ 
2004; ISommer-Larsen et al.ll2005l) . At early times the ICL mor- 
phology is dominated by long, linear features like streams that 
become more diffu se at later t imes as they spread in the clus- 
ter volume dRudick et al.l2009h . Most of the simulated intraclus- 
ter stars become unbound from their parent galaxies during the 
merging history leading to the formation of the brightest clus- 
ter galaxy (BCG) in the cluster center, while for ICL at larger 
radii other mechanisms like tidal stripping are more important 
dMurante et~aT1l2007l: IPuchwein et aTll2010h . 



In this paper we report surface photometry and long slit 
spectroscopy of the ICL in the center of the Hydra I cluster, a 
medium compact cluster at ~ 50 Mpc distance in the Southern 
hemisphere, with a central cD galaxy, NGC 3311. The aim is 
to compare the structures in the surface brightness distribution 
in the cluster core around NGC 3311 with kinematic informa- 
tion, including also the LOS velocities of intracluster planetary 
nebulas (ICPNs). Studying the kinematics of the ICL in nearby 
clusters like Hydra I is possible with ICPNs because (i) these 
objects are rel atively easy to detect due to their strong [01 1 1] 
emission line dJacobvll 19891 ICiardullo et al.ll 19891) . and (ii) they 
are goo d tracers of the light distribution of the par ent stellar pop- 
ulation dBuzzoni et alJl2006t ICoccato et alJl2009h . The LOS ve- 
locity distribution (LOSVD) of the ICPNs associated with the 
diffuse lig ht in the central 100 kpc arou nd NGC 331 1 was deter- 
mined by IVentimiglia et alJ (120081 1201 ll) . They found discrete 
velocity components at redshift ~ 1800 km s~' and ~ 5000 
kms -1 , in addition to a broad component with velocity dis- 
persion cr s 500 km s -1 approximately at th e systemic veloc- 
ity of the Hydra I cluster 3900 km s~'). Ventimiglia et alJ 
d201 lb concluded that the broad velocity component in the ICPN 
LOSVD m ay trace the high-veloci ty dispersion outer halo of 
NGC 3311 d Ventimiglia et alj|2010h . while the discrete compo- 
nents trace sub-components that fell through the cluster core, 
were tidally disrupted, and have not yet dynamically mixed in 
the gravitational potential. 



This paper is structured as follows: in Section [2] we present 
optical V-band images for the Hydra I cluster core. Isophote fit- 
ting and analysis of the surface brightness profiles is carried out 
in Section [3] Two-dimensional models for the optical data are 
derived in Section |4] showing the existence of an off-centered 
outer halo around NGC 3311 as well as tidal streams super- 
posed on this halo. In Section|5] long slit spectroscopic data and 
kinematic measurements for the halo of NGC 3311 and the lo- 
calized stream about HCC 026 are presented and discussed. In 
Section[6]we investigate the correspondence between the photo- 
metric components and kinematic substructures in the velocity 
distribution of Hydra I PNs, dwarf and SO galaxies within 50 kpc 
of the center of NGC 331 1. Section[7]briefly discusses the pecu- 
liar outer halo of NGC 331 1, the properties of the newly discov- 
ered tidal streams, and the formation of ICL from this group of 
galaxies in disruption. Finally, Section [8] contains our summary 
and conclusions. We assume a distance to the Hydra I cluster of 
D = 50 Mpc, so 1" = 0.247 kpc. 



2. Optical imaging of the Hydra I cluster core 

The core of the Hydra I cluster is dominated by two giant el- 
liptical galaxies, NGC 3311 and NGC 3309. Early CCD sur- 
face photometry showed that both NGC 3311 and NGC 3309 
are fitted by an R 1 ^ 4 la w within 30" distance from their centers 
dVasterberg et al.l 99 ll). On the bas i s of th e large R e = 98" value 
from the R 1/4 fit. IVasterbergetall (fl99lb classified NGC 3311 
as a cD galaxy. In the center of this galaxy, a complex dust 
lane is found, and the observed surface brightness profile thus 
falls below the best fit R 1 ' 4 profile. Both the dust lane and the 
core in the central pro file are confirmed by HST WFPC2 imag- 
ing dLaine et alJl2003h. NGC 3309 does not have a d ust lane in 
its core dVasterberg et al.lll991 ). Misgeld et al. (2008) published 
VLT/FORS 1 i mages of the Hydra cluster in the V-band, that 
were used by Rich tler et alJ (1201 ll) to derive the average radial 
surface brightness profile for NGC 3311 down to fiy = 26.5 mag 
arcsec -2 . They classified NGC 3311 as a cD galaxy because of 
the low central surface brightness and extended radial profile. 

In the present work we derive accurate quantitative photom- 
etry for the galaxy NGC 331 1, NGC 3309 and the diffuse light 
in the Hydra I core. In what follows, we use two data sets for 
this purpose: VLT/FORS1 archive data for the center of the 
Hydra I cluster and V band images obtained with WFI at the 
ESO/MPI 2.2m telescope. These V band images allow us to 
measure the profile shape at large radii where the surface bright- 
ness is low, and analyze the substructures in the light distribu- 
tion. The VLT/FORS1 and 2.2m/WFI data contain complemen- 
tary information: the FORS1 data are deep, and allow reliable 
measu rements down to fiy = 26.5 mag arcsec -2 dRichtler et al.l 
120111) . while the WFI data with their large field-of-view provide 
a robust estimate of the sky background. 



2. 1. VLT/FORS1 deep V band photometry - Observations 
and data reduction 

Johnson V band imaging data for NGC 3311 and NGC 3309 
were retrieved from the ESO Science Archive Facility, acquired 
with VLT/FORS 1 in service mode on April 04, 2000 (program 
65.N-0459(A), PI Hilker). The VLT/FORS 1 images cover a 
field-of-view (FoV) of 6.8' x 6.8' at an angular scale of 0".2 
pixel -1 . During the observations, the FoV was centered at a = 
10h36m36.14s, 6 = -27d32m51s, avoiding a bright star NE of 
NGC 3311. The dataset included three exposures of 480 sec each 
for a total observing time of 0.4 hr. The resulting average seeing 
in the combined median image is FWHM ~ 0".6; the central 
~ 5" diameter region of NGC 3311 and ~ 8" diameter region of 
NGC 3309 are saturated. 

Standard calibrations, bias, and sky flats were also retrieved 
from the ESO Archive. Several Landolt standard stars in the 
Rubin 152 field were observed in V band for the photometric 
calibration. The zero point for the V band VLT/FORS 1 pho- 
tometry ZPy fORsx = 27.43 + 0.06 is derived independently in 
the present work. The correction for extin ction in the V band 
amounts to 0.25 mags dSchlegel et al.lll998l) . 

The instrumental signature of the VLT/FORS 1 images was 
removed using standard IRAF tasks. Spikes from saturated stars 
in the field were also removed; see Section 12.21 Images were 
registered before co-addition; the co-added final VLT/FORS 1 V 
band image is shown in Fig.Q] 



2 



M. Arnaboldi et al.: Two tails and an off-centered envelope in the Hydra I cluster core 



* - 



50" 



• 


■ 

■ 

* 

V 


14 
1 

.. ! 


* 














50' 











• " • N 






• 


* 






SO" 







Fig. 1. Left panel - VLT/FORS1 V band co-added image of the Hydra I cluster core with NGC 331 1 and NGC 3309 visible in 
the upper part of the image. A 2D-photometric model is fit within the region of 6'. 8 x 6' .4 limited by the box (see Sect. 14. II for 
details). Middle panel - The mask "cmask" adopted for the 2D-model fit of the VLT/FORS 1 data, which masks stars and galaxies in 
the chosen region as well as the centers of NGC 331 1 and NGC 3309. Right panel - The mask "allmask" for the 2D fit that masks 
also part of the North-East regions in addition to the stars and galaxies in "cmask" (see Sect l4.1l for details). In the middle and right 
panels, the two crosses indicate the positions of the centers of the two giant galaxies. 



2.2. Wide Field Imager Vband photometry - Observations 
and data reduction 

Johnson V band imaging of the Hydra I cluster was acquired 
in service mode on the night of January 12, 2006 at the Wide 
Field Imager (WFI) on the ESO/MPI 2.2m telescope, at the La 
Silla observatory. The WFI imager covers a field-of-view (Fo V) 
of 34' x 33' with a mosaic of 4x2 CCDs (2kx4k) at an angular 
scale of 0".238 pixel -1 . During the observations, the FoV was 
centered at a = 10h36m51s, 6 = -27d31m35s. 13 exposures of 
300 sec each were obtained for a total observing time of ~ 1.1 
hr. The resulting average seeing in the combined median image 
is FWHM ~ 0".7. 

Standard calibrations, bias, sky flats and dark skies, were 
also obtained. Several Landolt standard stars in the Rubin 149 
field were observed in V band for the photometric calibration. 
The zero point for the V band photometry is ZPy = 24.02+0.02. 

Data reduction is carried out with standard IRAF tasks for 
pre-reduction and calibration. After bias subtraction and flat 
fielding, the frames are corrected for any residuals by using dark 
skies. Then the average background emission is measured in sev- 
eral regions of the FoV far from the galaxy light and the final av- 
erage value is subtracted off each single frame. The IRAF task 
NOAO.NPROTO.IRMOSAIC is used to obtain a mosaic of the 
CCDs. Each frame is inspected and bright spikes associated with 
saturated stars are removed by linear interpolation from nearby 
columns/rows before frames are co-added. Residual sky offsets 
with respect to the average value are removed before the CCD 
frames are mosaiced. Finally, the mosaic images are registered. 
Then image fluxes are scaled to a reference image to account for 
transparency variations, following which images are co-added 
with a threshold rejection implemented to reduce effects from 
CCD gaps. The final co-added image is shown in the left panel 
of Fig. 

As a first step in the study of the light distribution in 
NGC 3311, we determine the extension of the dust lane in 
its central region. We use the FMEDIAN task in IRAF with a 
smoothing box of 15 x 15 pixels, and compute the ratio of the V 
band image and its FMEDIAN smoothed version; the resulting 
unsharp masked V-band image is shown in Figure [3] 

This figure illustrates the presence of several distinct compo- 
nents at the center of NGC 3311. A complex dust lane crosses 



the galaxy center in the direc tion North-Sou th (NS) in the high 
angular resolution image in lLaine et al.l (120031) . Bright knots are 
seen East and Southwest of the galaxy center, wit hin and around 
the du st lane out to 8" in radius. According to Vasterberg et al.l 
d!991h . the dominant knot is bluer, A(B — r) = -0.10, than the 
surrounding stellar population. 




Fig. 3. Unsharp masked V band image obtained from optical data 
acquired at the ESO/MPI 2.2m telescope with the WFI; North is 
up, and East to the left. The white bar indicates 10" length. The 
inner dus t lane at the center of NGC 331 1 is about 2" wide and 
7" long (lLaine et al 1 12003b . and is embedded within a central 
light excess of about 16" in diamet er, where bright k nots are 
visible here and in the HST image of lLaine e t al. (2003). 



3. Surface photometry for NGC 3311 in V band 

In this Section, we describe isophotal analysis for the light dis- 
tribution of NGC 3311 from the V-band data. We then derive 
V-band surface brightness profiles along several specific direc- 
tions through the galaxy to characterize the deviations from a 
simple concentric elliptical light distribution. 
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Fig. 2. Left panel - 2.2m/WFI V band image of the Hydra I cluster core; NGC 33 1 1 is visible on the left and NGC 3309 on the right. 
A 2D-photometric model is fit within the region of 6'. 4 x 4'. 8 limited by the box. Right panel - The mask "allmask" adopted for the 
fit to the bright central regions (see Sect. 14.2 I for details). The two crosses in the right panel indicate the positions of the centers of 
NGC 3311 and NGC 3309. 



3. 1. Isophote fitting 

We use the ELLIPSE task in IRAF on the V images to perform 
the isophotal analysis of NGC 331 1. Position angle (P.A.), ellip- 
ticity, and average surface brightness profiles for NGC 331 1 are 
shown in Figure [4] 
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Fig. 4. One dimensional profiles of the RA. for the isophote 
major axis, ellipticity, and surface brightness obtained with 
ELLIPSE from V-band images of NGC 331 1 (VLT/FORS1 blue 
crosses, 2.2m/WFI black open dots). Because of the central dust 
lane, the profiles are well defined only outside 3 ".5. 



The ELLIPSE V-band average surface brightness profile ex- 
tends out to 80" from the galaxy center before the halo light from 
NGC 3309 becomes significant. Note the excellent agreement 
between the VLT/FORS1 and WFI surface brightness profiles. 



The total integrated magnitudes inside circular apertures, whose 
radius corresponds to the last measured point, are: niy FI (R = 



77" .4) = 11.55 mag; m 



FORS 



71".5) = 11.36 mag. The 



the last measured point are Rjy' 
r fors _ 32 » + 2 



half-light radii evaluated from the light-growth curves out to 

= 36" + 2" (~ 8.9 kpc), 

(~7.9kpc)." 

At R < 5", the presence of a dust-lane is signaled by the large 
variations in the V-band ellipticity and R A. profiles. Between 5" 
and 15" the V-band isophotes twist by about 20 degrees: this is 
caused by the presence of the inner bright knots seen around the 
dust lane (see Sect. |2,2J . A twist of the NGC 3311 isophotes is 
also reported by Vasterberg et al. (1991), but the RA. variation 
was not quantified. For radii 15" < R < 45", the V-band ellip- 
ticity and RA. profiles are nearly constant, with average values 
e - 0.05 and RA.^ 32°, i.e., the isophotes in this region are 
nearly round. 



3.2. Analysis of surface brightness profiles along principal 
axes of the light distribution 

The position angle profile from the isophote fit in Fig. |4] has 
a complicated radial dependence which we show now is re- 
lated to asymmetries in the light distribution. To this end we 
investigate separately the surface brightness profiles along sev- 
eral axes across the light distribution in NGC 3311. From the 
isophote fitting, the major axis is at P.A.= 32° and the minor 
axis is at P.A.= 122°. Along these directions, we have extracted 
one-dimensional (ID) surface brightness profiles from V-band 
VLTVFORS 1 and 2.2m/WFI images. 

The V band surface brightness profiles become asymmet- 
ric for R > 30", as shown in Figs. [5] and [6] The major-axis 
profile at P.A.= 32° (North East - NE) is up to one magnitude 
brighter than that at RA.= 212° (South-West - SW) in the radial 
range 20" < R < 120" (Fig. [5]). There is additional light along 
P.A.= 212° at R > 100" which may be related with a substruc- 
ture in the Hydra I core. This will be further investigated in the 
residual image, once the light distribution from the bright galax- 
ies is accounted for. The light profiles along the minor axis, at 
P.A.= 122° (South East - SE) and P.A.= 302° (North West - NW) 
are symmetric for R < 60" (Fig.©. At RA .= 302° and R > 60" 
we see the light from the NGC 3309 halo. For R < 5", both V 
band profiles show the absorption by the central dust-lane. 
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Fig. 5. Upper panel: folded V-band profiles extracted along 
P.A.= 32° (major axis, 2.2m/WFI black full lines, VLT/FORS1 
blue dashed lines). NE is along P.A.= 32° and SW along P.A.= 
212°. Lower panel: difference profiles. The folded profile along 
the major axis illustrates the excess of light in the NE quadrant 
of the NGC 3311 halo, in the range of radii 20" < R < 120", 
with a maximum excess of about one magnitude at R = 80". 
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Fig. 6. Same as Fig|5] Upper panel: folded V-band profiles ex- 
tracted along P. A.= 122° (minor axis, 2.2m/WFI black full lines, 
VLT/FORS1 blue dashed lines). SE is along P.A.= 122°, NW 
along P.A.= 302°. Lower panel: difference profiles. On the NW 
side, for R > 65" we see the contribution from the outer regions 
of NGC 3309. 



In addition, we extract folded profiles along the NS (P.A.= 
0°) and EW (P.A.= 90°) directions. These are also asymmetric; 
see Fig. □ and Fig. [8] At P.A.= 0° and R > 20", the galaxy 



Fig. 7. Same as Fig[5] Upper panel: folded V-band profiles ex- 
tracted along P.A.= 0° (2.2m/WFI black full lines, VLT/FORS1 
blue dashed lines). N is along P.A.= 0°, S along P.A.= 180°. 
Lower panel: difference profiles. 
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Fig. 8. Same as Fig|5] Upper panel: folded V-band profiles ex- 
tracted along P.A.= 90° (2.2m/WFI black full lines, VLT/FORS 1 
blue dashed lines). E is along P.A.= 90°, W along P.A.= 270°. 
Lower panel: difference profiles. 



becomes brighter than at P.A.= 180°; the difference is about 0.5 
mag at R = 80" (Fig. 01. At P.A.= 90°, 270° the profiles are 
symmetric for R < 30", while in the radial range 30" < R < 60" 
the profile at P.A.= 90° is brighter (by about 0.2 mag at R ~ 50") 
than that at P.A.= 270°. At larger radii the profile is affected by 
the light of NGC 3309. Figs.Qand [8]illustrate that the additional 
light is distributed over an opening angle larger than 90° as seen 
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from the center of NGC 3311, therefore the substructure leading 
to the asymmetry in Fig.|5]is not a radial stream or fan. 

The VLT/FORS1 and 2.2m/WFI profiles in Figs. and 
show very good agreement on both the N/NE and S/SW sides 
of NGC 3311 and also in the difference profiles. This shows 
that systematic effects caused by scattered light from the bright 
star NE of NGC 331 1 are negligible and that there are no resid- 
ual gradients due to the sky background along these directions. 
Figs.[6]and[8]show a small relative variation between the W/NW 
and E/SE sides, amounting to ~ 0.2 mag over 200". As we will 
discuss further in Sect. 14.21 this does not affect our main results. 

The analysis presented in this Section illustrates that the halo 
of NGC 3311 is asymmetric. Between P.A. = 0° and P.A.= 90°, 
and in the radial range 30" < R < 120", it is brighter than in the 
remaining three quadrants. 

4. Models for the light distribution in NGC 331 1 , 
NGC 3309, and the residual diffuse light in the 
Hydra I cluster core 

In this Section we construct two dimensional (2D) surface 
brightness models for the light distribution in the Hydra cluster 
core around NGC 331 1. Any such model must include the light 
distribution of NGC 3309, because the outer regions of the two 
elliptical galaxies overlap along the LOS. Furthermore, the mod- 
eling must also account for the asymmetric bright component in 
the NE outer region of NGC 3311, signaled by the asymmetry 
of the extracted profiles in Section l3~2l 

4. 1. V-band 2D model for NGC 331 1 and NGC 3309 

We model the surface brightness distribution of NGC 331 1 and 
NGC 3309 from the V band VLT/FORS1 image. For fitting the 
model to the VLT/FORS 1 data, we first adopt the mask shown in 
the middle panel of Fig.[TJ identified as "cmask" in what follows. 
It masks the saturated/bright stars in the field, the central R < 5" 
of NGC 3311, a region R < 8" around the saturated center of 
NGC 3309, and the background galaxies. 

For the "cmask" model fit, the 2D Sersic model has param- 
eters n v = 3.3, R e = 21". 1 for NGC 3309 and n v = 5.0, 
R e = 270".6 for NGC 3311; see Tables Q] and The residual 
V band image shows a number of structures. In particular, to the 
SW of NGC 3311 around the galaxy's major axis at P.A.= 212°, 
in a ~ 90°-wide cone and radial range 40" < R < 80", we 
see (unphysical) negative residuals of about 260 ADU per pixel 
with respect to the sky counts measured in the "empty" eastern 
boundary of the VLT/FORS 1 image. This negative residual sig- 
nals that NGC 3311 is fainter than the 2D Sersic model there. 
It derives from fitting an azimuthally averaged model to an in- 
trinsically asymmetric light distribution; we saw previously in 
Sect. 13.21 that the galaxy is brighter in the NE quadrant than in 
the other three quadrants. 

Therefore, we next investigate whether by masking the 
whole North-East quadrant of NGC 3311, the "symmetric" 2D 
model for the V band light in NGC 3311 will acquire a lower 
value for the Sersic index «, and thus the problem of negative 
residuals at P.A.= 212° is reduced. We define a new mask which 
additionally covers most of the NE quadrant in the NGC 3311 
halo: this is shown in the right panel of Fig. [JJand is identified as 
"allmask". When the light excess is masked, the GALFIT 2D fit 
to the V band FORS1 image returns a Sersic index «y = 4.8 and 
R e = 198".8 for NGC 3311. The GALFIT model of NGC 3309 
and NGC 331 1, using "allmask" with the VLT/FORS 1 data, and 



Parameter V band V band V band 

FORS1 FORS1 FORS1 

max. symmetr. 



Comp. type 


Sersic 


Sersic 


Sersic 


m tot 


11.6 ±0.05 


11.6 ±0.05 


11.2 ±0.01 


R, 


21". 1 ±0".07 


21".9±0.07 


35" ± 0.25" 


n 


3.3 ± 0.02 


3.3 ± 0.02 


5.7 ±0.04 


b/a 


0.87 ± 0.01 


0.85 ± 0.01 


0.87 ± 0.01 


P.A. 


48.2 ± 0.2 


50.7 ± 0.2 


49.8 ± 0.2 


mask 


"cmask" 


"allmask" 


"allmask" 



Table 1. Parameters for the 2D GALFIT models of NGC 3309. 
Apparent V band total magnitudes are in agreement with the 
published NED values (RC3). 



the residual image (difference between image and model) are 
shown in Figure [9] Parameters are given in Tables [T]and|2] 

The residual image still shows negative values in the same 
region in a ~ 90°-wide cone around P.A.= 212° and for the ra- 
dial range 40" < R < 80", but the level is now reduced to 110 
ADU per pixel, compared to 260 ADU per pixel in the "cmask" 
residual image. Therefore the extra light in the NE of NGC 3311 
must extend slightly further than the masked area. 

In addition, the residual image shows several positive struc- 
tures which correspond to additional light in the real image with 
respect to the 2D models of NGC 33 1 1 and NGC 3309, at a level 
of 200-500 ADU per pixel relative to sky. These structures will 
be more quantitatively discussed in Sects. 1431 and !4.3l in the con- 
text of a "maximal" symmetric model with which the negative 
residuals are removed. Briefly, 

- In the outer regions, NGC 3311 is more luminous than the 
best fit 2D model in a wide area covering the whole NE quad- 
rant to the galaxy. Superposed on this envelope we also see 
a prominent tidal stream around HCC 026, one of the dwarf 
galaxies in the Hydra I cluster core. The stream is outlined 
by the red dashed lines in Fig.|9]around ~ 43" distance from 
the NGC 3311 center. The tidal tails around HCC 026 extend 
about 15 kpc NW and 5 kpc NE of the dwarf galaxy. 

- The SO galaxy HCC 007 situated 1 19" South of NGC 331 1 
appears to be embedded in an extended, thick, tidal stream 
visible at faint surface brightness levels. Emerging from 
HCC 007, the stream circles the south part of the Hydra I 
core, it encompasses both NGC 3311 and NGC 3309, and 
extends over at least 450" (~ 110 kpc), as indicated by the 
dashed blue lines in Fig. [9] The eastern tail of HCC 007 ap- 
pears to join the bright envelope to the NE of NGC 33 1 1 ; as 
part of the envelope it may continue considerably further to 
the north. 

Fig.flOlshows ID surface brightness profiles extracted along 
the NE and SW sides of the major axis of NGC 3311, as illus- 
trated in Fig. [9] For radii 5" < R < 30", the extracted profiles are 
straight lines in the R 1 ^ 4 plot: deviations from a de Vaucoleurs 
profile become large at R > 30" on the NE side (P.A.= 32°), 
indicating additional light there. The flattening of the surface 
brightness profile at P.A.= 32° in Fig. [10] with respect to the 
steeper profile at P.A.= 212° indicates that at radii larger than 
30", the outer surface brightness contours shift towards the NE; 
i.e., the envelope is off-centered to the NE relative to the bright 
central region of NGC 3311. This confirms the asymmetries of 
the extracted surface brightness profiles discussed in Sect. [3] 

Fig. [10] also explains the origin of the negative residuals 
in the "allmask" model: the surface brightness profile fitted by 
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10' 






Fig. 9. GALFIT model of the VLT/FORS 1 image obtained with "allmask" (see right panel of Fig. Left panel: VLT/FORS 1 V 
band image. Full and dashed red lines indicate the major axis at PA=32°, and the minor axis, respectively. Center panel: combined 
2D GALFIT model for NGC 3311 and NGC 3309. Right panel: residual image, showing extra light associated with the NE outer 
halo of NGC 331 1, a region of negative residuals ~ 60" South West of NGC 331 1, and preliminary evidence for two tidal streams 
around HCC 026 and HCC 007, outlined by the red- and blue dashed curves. The image size in all panels is 6'. 8 x 6'. 4, and darker 
colors indicate brighter regions. 



Parameter 


Vband 


Vband 


Vband 




FORS1 


FORS1 


FORS1 








max. symmetr. 


Comp. type 


Sersic 


Sersic 


Sersic 


m ',o, 


9.7 ± 0.01 


10.1 ±0.01 


9.52 


K 


270". 6 ± 3" .4 


198".8±2".2 


850" ± 50" 


n 


5.0 ± 0.02 


4.8 ± 0.02 


10.5 ±0.1 


b/a 


0.93 


0.93 


0.93 


P.A. 


38 ± 0.5 


32 


32 


mask 


"cmask" 


"allmask" 


"allmask" 



Table 2. Parameters for the 2D GALFIT fit models of 
NGC 3311. 



* The nitot and R e values computed from the GALFIT fit model 
for NGC 3311 are respectively brighter and larger than those 
computed in Sect. 13. II from the V band light growth curve, be- 
cause the latter is measured only out to 71".5, while the 2D 
model is fit to the whole 6'. 8 x 6'.4 area. 



GALFIT (green line) falls slightly above the SW (P.A.= 212°) 
major axis profile because of the surface brightness drop around 
30 arcsec, so that the symmetric 2D model built from this pro- 
file generates negative values in the residual image about P.A.= 
212°. 

In the next Section, we will improve on this by construct- 
ing a "maximal" symmetric model for the light centered on 
NGC 331 1. Then we will use this to construct a new 2D model, 
and to discuss the structures in the corresponding residual image 
in more detail. 

4.2. The asymmetric light distribution around NGC 3311: 
extended off-centered envelope and tidal tails 

We build a maximal symmetric model (MSM) for NGC 3311 
in two steps. First, we fit the steep surface brightness profile to- 
wards SW, along P.A.= 212°, with a Sersic law. Second, we build 
a 2D GALFIT model with the Sersic parameters fixed by the ID 
fit to the steep profile, which therefore leads to positive or null 
residuals everywhere in the NGC 331 1 halo. 

Fig. [10] shows the ID fit profile over-plotted on the data 
(red line); the model fits the data well in the radial range 15"- 
80". The parameters are n v = 10.5 + 0.1, R e = 850" + 50", 




R 1 ^ 4 (arcsec) 



Fig. 10. Folded ID surface brightness profile along the major 
axis of NGC 3311 plotted as function of R 1 ^ 4 . Open symbols 
show values along P.A.= 212° (SW), full symbols those along 
P.A.= 32° (NE of the galaxy center). The green continuous line 
shows the extracted ID major axis profile of the Sersic "allmask" 
model for the V-band FORS1 data described in Sect. 07TJ The 
red continuum line is the major axis profile for the light centered 
on NGC 3311 ("maximal" symmetric model, MSM). The blue 
full line shows the ID fit to the surface brightness profile along 
P.A.= 32° (NE) with a two component model, consisting of the 
MSM plus a one-sided exponential profile . The arrows indicate 
the position of the stream around HCC 026 and the light at R 
80" associated with the dwarf galaxy HCC 027 which falls on 
the extracted profile. See Sect. l4.2l for details. 
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and yU e = 28.07 + 0.1 mag arcsec" 2 ; see Table [2] Now we 
use GALFIT with the mask "allmask" to generate a 2D MSM 
for NGC 3311, using the Sersic law with parameters set from 
the ID best fit, and the major axis PA=32° and axis ratio 0.93 
from our previous 2D model fits. Then, keeping this component 
fixed, we use GALFIT again to determine new best fit param- 
eters for NGC 3309, matching a Sersic model to the remaining 
light. Results of this NGC 3309 fit are listed in Table Q] Adding 
both components, we obtain the MSM for the two bright galaxies 
NGC 331 1 and NGC 3309. By construction, the residual image 
obtained as the difference between the V-band VLT7FORS 1 im- 
age and this MSM should show only positive or null residual 
values around NGC 3311 (outside the central R < 15" contain- 
ing dust and bright knots). 

This residual image, displayed in Fig. Q~T] clearly shows the 
presence of an additional, extended envelope, off-centered to- 
wards the NE of the inner regions of NGC 3311. It is high- 
lighted by the white dashed ellipse in Fig. QT| The residual im- 
age also shows the bright tail around HCC 026, and the broader 
tail around HCC 007 which connects with the off-centered en- 
velope. These structures stand out more clearly in Fig. HTIthan 
in the earlier Fig. |9]because with the MSM we no longer have 
negative residuals. Note that the "half -ring" morphology around 
NGC 331 1 slightly inside the HCC 026 stream probably signals 
that a small fraction of the light in the "maximal" symmetric in- 
ner halo for r < 30" should be counted as part of the off-centered 
envelope, such that its surface brightness on top of NGC 3311 
becomes a constant /iy = 23.4 mag arcsec -2 . However, we have 
no independent way to determine this light fraction. The "maxi- 
mal" symmetric model is constructed such as to put the minimal 
fraction of light into the off-centered envelope. 

To quantify the asymmetry in the NGC 3311 halo, we also 
fit the NE surface brightness profile in Fig.[l0](P.A.= 32°) with a 
two-component model, consisting of a Sersic law with the same 
n and R e and an additional exponential profile for the outer re- 
gions. The exponential profile often is a goo d approximation fo r 
the outer stellar envelopes of cD galaxies (Sei gar et al.l 12007). 
This combined model provides a good fit to the NE profile 
which is also shown in Fig. [10] (blue line); the parameters are 
n v = 10.5 + 0.1, R e = 850" + 50" and fi e = 28.2 + 0.1 mag 
arcsec -2 for the Sersic component, and fio = 23.2 + 0.1 mag 
arcsec -2 and r;, = 200" + 50" for the exponential component. 
The peak surface brightness of the off-centered envelope in the 
"maximal" symmetric halo model is //o,v = 23.2 mag arcsec -2 , 
as given by the exponential fit along P.A.= 32°. 

In the residual map, Fig.[TT| the measured surface brightness 
at the position of the HCC 026 tail is also ~ 23.2 mag arcsec -2 . 
This is the combined surface brightness of the off-centered enve- 
lope and the tail, which are superposed along this line-of-sight. 
The contribution of the HCC026 tail within the polygon area 
limited by the blue dashed line in Fig. QT|is about 23%; see 
the next section. Hence the average surface brightness of the tail 
around HCC 026 in the area enclosed by the blue dashed poly- 
gon in Fig.Q~T]is ftv = 24.8 + 0.2 mag arcsec -2 . 

The surface brightness of the tail around HCC 007 can be 
directly measured on the residual image. It is brighter on the NE 
side where fiy = 24.4 ± 0.5 mag arcsec -2 , than on the NW side 
where it is about one magnitude fainter. The brighter NE tail 
of HCC 007 is clearly a localized structure with higher surface 
brightness than other areas at similar distance from either NGC 
33 1 1 or NGC 3309. The central region of the brighter NE tail of 
HCC 007 contains two small objects that appear to be connected 
by a possible further small stream. 



We can check with Figure [12] that the NW tail of HCC 007 
is not an artifact of residual halo light in NGC 3309. The up- 
per panel of this figure shows that the minor axis profile of this 
galaxy (P.A.= 320°) is well described by the fitted Sersic pro- 
file, whereas additional light is seen along P.A.= 230° and par- 
ticularly along P.A.= 190° which goes through the brightest part 
of the NW tail close to HCC 007. However, it cannot be ruled 
out that this region contains some light from a very diffuse, low 
surface brightness structure around NGC 3311 and the cluster 
center. 

Finally, we show that the morphologies of the substructures 
seen in the residual image in Fig. QT|are n °t due to systematic 
effects in the 2D modeling caused by uncertainties in the back- 
ground estimates. This concern can be addressed by showing 
that these substructures are present in two independent datasets, 
the VLT/FORS1 and 2.2m/WFI images. We have already seen 
in Sect. I3.2l that the asymmetries in the extracted profiles along 
P.A.= 32°, 212° are very similar in the two data sets; this con- 
firms independently the presence of the off-centered envelope 
NE of NGC 3311. The lower panel of Fig. [12] compares the 
VLT/FORS1 and WFI surface brightness profiles through the 
fainter NW tail of HCC 007 along P.A.= 190° as seen from the 
center of NGC 3309. The light excess in the HCC 007 tail is clear 
in both data sets. Both results are independent of the slight resid- 
ual gradient described in Sect. 13. 21 because this is essentially per- 
pendicular to P.A.= 32° and P.A.= 190°. Finally, we have also 
carried out 2D GALFIT modeling on the 2.2m/WFI V band im- 
age, using the ID maximal model fit parameters for NGC 331 1, 
and the "allmask" mask shown in Fig.[2]to determine the param- 
eters of NGC 3309. The residual V-band 2.2m/WFI image has a 
lower S/N than that based on the VLT/FORS1 data. Nonetheless 
both elongated streams around HCC 026 and HCC 007 are inde- 
pendently confirmed. In the case of HCC 007, only the brighter 
NE tail is visible in the residual image. 

4.3. Luminosities of the substructures in the Hydra I core 

We next compute the total luminosity of the off-centered enve- 
lope from the flux within the elliptical aperture shown in Fig.fTTI 
We mask stars, dwarf galaxies, and a circular region with R = 
25" centred on NGC 3311. The resulting luminosity in the off- 
centred envelope is Lv,NE,env = 1 -2 x 10 10 (+6.0 x 10 s ) LgQ using 
a distance of 50 Mpc for the Hydra I cluster. This estimate does 
not include the luminosity of the MSM halo in this region, but it 
does include the light in the stream around HCC 026, which is 
seen in addition to the exponential component in Fig.fTOl 

For comparison, we compute the luminosity of the MSM 
i) within the same elliptical aperture, and ii) in the circular 
annulus between R = 25" and 120", by integrating the flux 
of the GALFIT model output image over the respective areas. 
The luminosities of the MSM model evaluated in these two re- 
gions are Lv.NE.haio = 2.3 x 10 10 (+6.0 x 1O 8 )L and Lv.haio = 
7.5 x 10 10 (±2.6 x 10 9 )L o . Thus, within the area limited by the 
elliptical aperture in Fig. QT| the luminosity in the off-centered 
envelope amounts to 50% of the light in the MSM, while it is 
only 15% of the total luminosity of the MSM integrated between 
25" and 120". 

To estimate the luminosity of the HCC 026 tail we carry 
out surface photometry in a polygon designed to include most 
of the NW half of the stream; it is constructed with the IRAF 
POLYMARK and POLYPHOT tasks and is shown by the blue 



1 The 1 - tr error in the luminosity is estimated by integrating the sky 
rms per pixel over the area enclosed by the elliptical aperture. 
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26.93 25.26 24.63 24.24 23.95 23.72 23.53 23.37 23.23 

Fig. 11. Residual VLT/FORS1 V-band image with respect to the MSM for NGC 331 1 and NGC 3309. The white ellipse indicates 
the extended halo off-centered towards the NE of the bright inner regions of NGC 3311 (black). The angular extent of this envelope 
as seen from the center of NGC 331 1 is consistent with the opening angle of the additional light as shown by Figs. Q and [§]. The 
most noticeable substructure on top of the extended envelope NE of NGC 331 1 is the tidal stream associated with the dwarf galaxy 
HCC 026. Its NW part is indicated on the image with the blue dashed polygon. At larger radii, and at fainter surface brightness 
levels, we see a thicker tail around HCC 007; the luminosity in this tail is computed in the areas limited by the white polygons, see 
Sect. 14.31 The two small objects near the center of the NE polygon appear to be connected by a further small stream; see also Fig. [9] 
The FoV is 6'. 8 x 6'.4. Colour levels are labeled with the corresponding surface brightness values. This image shows that in the 
region south of HCC 007 the background is flat, with a systematic E-W gradient across the image at less than ~ 10 ADU. The white, 
green, and blue lines from the center of NGC 3309 indicate the P.A.s at = 190°, 230° and 320° respectively; see extracted profiles in 
Fig. 03 



dashed polygon in Fig. QT| The polygon defined for the 
HCC 026 tail is at an average distance of 53". From the lu- 
minosity measured in this polygon we subtract the contribu- 
tion of the smooth envelope in the same area, estimated from 
the average residual surface brightness on both sides of the 
stream. This gives us a differential measurement of the lumi- 
nosity of the HCC 026 tail in the polygon. The resulting lu- 
minosity of the NW part of the stream is Lv.NWstreamHCC 026 = 
4.8 x 10 8 (+8 x 10 7 )L o , which is equivalent to an apparent mag- 
nitude of mv jS treamHCC 026 = 16.6. We find that the NW stream 
of HCC 026 contributes about ~ 15% of the combined sur- 
face brightness of the MSM halo and off-centered envelope at 
the stream position (see also Fig. [Tot , consistent with the esti- 



mate obtained independently from stellar population analysis in 
ICoccato et al.l d201 lab . The luminosity in the southern part of the 
HCC 026 stream is more difficult to disentangle from the diffuse 
halo, but comparable. It is seen in Fig.[8]as an excess of surface 
brightness of ~ 0.1 mag over ~ 15" in the E profile over that on 
the W side. The total stream luminosity is clearly several times 
larger than the total luminosity of the dwarf galaxy H CC 026 it- 
self, L v ,hcc 026 = 1.5 x 10 8 L o dMisgeld et al.ll2008l) . and if the 
stream is physically related to HCC 026, a large part of the lu- 
minosity of this galaxy has already been tidally dissolved! 

In a similar way, we compute the light in the tail emerging 
from HCC 007. We define two polygons shown by the white 
dashed boxes in Fig.[TT]NE and NW of HCC 007 and carry out 
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Fig. 12. Upper panel: surface brightness profiles extracted in 
5° wide cones along the major (P.A.= 230°, green squares), 
minor axis (P.A.= 320°, blue triangles) and intermediate axis 
(P.A.= 190°, black line) of NGC 3309, plotted as function of 
R ll \ The PA. at 190°, 230° and 320° are indicated by the white, 
green, and blue lines from the center of NGC 3309 in Fig. Q~T] 
The red continuous line shows the extracted ID major axis pro- 
file for the 2D GALFIT "maximal" symmetric model (MSM) to 
the light of NGC 3309. While this is a good fit for the minor axis 
profile at large radii, profiles along P.A.= 230° and P.A.= 190° 
show additional light, as they go trough the brightest part of the 
HCC 007 NW stream. Lower panel: surface brightness profile 
from VLT/FORS1 (black line) and 2.2AVFI (red line) V band 
images along P.A.= 190°. 



the photometry with the POLYP HOT task. The polygons defined 
for the HCC 007 tails are at average distance of 120" and 173", 
where the surface brightnesses in the stream are 24.4 and 25.4 
mag arcsec -2 , and those of the "maximal" Sersic model are 24.3 
and 24.8 mag arcsec -2 . The apparent V band magnitudes in 
the two polygons are m v = 14.9 and 16.3 for the NE, NW 
sides, respectively. At a distance of 50 Mpc, the combined lu- 
minosity in the two polygons on the HCC 007 tails is there- 
fore Lv.streamHCC 007 = 2.95 x 10 9 (+5 x 1O 8 )L . For compari- 
son, the V band luminosity of HCC 07 from the Vo magni- 
tude measured in Mis g eld et all (|2008), V Q = 14.18 + 0.01, is 
Lv,hcc 007 = 4.5 x 10 9 L Q . Thus the ratio between the lumi- 
nosity in the tail polygons and of the SO galaxy HCC 007 is 
iv.streamHCC 007 / L H cc 007 = 0.66. This is obviously a lower limit; 
especially if the NE tail extends into the off-centered halo of 
NGC 331 1, its luminosity may be substantially larger than esti- 
mated with the NE polygon. Thus at least 40% of the original 
luminosity of HCC 007 have already been dissolved from the 
galaxy. 

In summary, the extensive photometric analysis of our V 
band images provides strong quantitative evidence for an ex- 
tended envelope around NGC 3311 off-centered towards the NE, 
and for luminous streams around HCC 026 and HCC 007. We 



now wish to establish the kinematic association of these diffuse 
structures to their galaxies, with the goal of understanding their 
origin. 

5. Long slit spectroscopy for the NGC 331 1 outer 
halo and for the tidal stream of HCC 026 

First we investigate whether the presence of the off-centered, ex- 
tended envelope around NGC 3311 is associated with any clear 
signatures in the kinematics obtained from deep long slit spectra. 
In Section IBTTI we discuss the mean LOS velocity and velocity 
dispersion profiles with respect to the center of NGC 331 1. We 
then describe the link between kinematic and surface brightness 
features in order to constrain the dynamics of the inner galaxy 
and the outer halo. In Section 15.21 we seek spectroscopic evi- 
dence for the presence of the tidal tails emerging from HCC 026. 

5. 1 . Kinematic signatures of the off-centered envelope 
around NGC 3311 

In IVentimiglia et alj d2010h we already showed that the halo 
around NGC 3311 for R > 15" is dynamically hot. The kine- 
matic data for NGC 3311 show an extremely rapid rise from 
galaxy-typical values of <x L osOK < 15") 170 km s -1 to veloc- 
ity dispersions cr LOS (50") > 400 km s -1 dominated by the clus- 
ter potential. This study was based on Gemini long slit spectra 
along P.A.= 63°, and deep VLT/FORS2 long slit spectra cen- 
tered on the dwarf galaxy HCC 026 with PA.=142°. More re- 
cently, the VL T7FORS1 spectra along P.A.= 29° published by 
iRichtler et al.l (1201 ll) confirmed the large outer velocity disper- 
sion values and the asymmetry of the Vlos profiles, but also 
showed deviations between measurements at different PA. 

The mean LOS velocity VlosC^) profiles along both PA.s 
increase outwards from NGC 3311 by ~ 100 km s -1 to the 
NE, while they remain close to the galaxy systemic velocity 
VWy nW = 3800 km s -1 to the SW. Thus IVentimiglia et al.l 
(1201 Ol) concluded that the NE outer halo has slightly more red- 
shifted mean LOS velocities, VWo - 3900 km s , and is offset 
by - +100 km s -1 relative to the NGC 3311 center. The photo- 
metric evidence for an offset halo component to the NE of NGC 
331 1 shown in Fig. [TT] suggests that the stars in the offset halo 
could be responsible for the shift in LOS velocities. 

Building on the photometric model derived in Section 14.21 
we can check whether either (Tlos(R) and/or Vlos(R) increase 
at positions where the surface brightness of the extended offset 
envelope is comparable to or brighter than the surface brightness 
of the MSM for NGC 331 1. In Fig. [13] we illustrate this analy- 
sis, comparing kinematic and surface brightness profiles along 
the galaxy's approximate major axis (P.A.= 29°, and P.A.= 32°, 
respectively). Indeed, both Vlos and <x LO s increase as the sur- 
face brightness of the offset envelope increases from R st 15" 
outwards. In particular, the Vlos values are more red-shifted on 
the NE side of the galaxy center than either at the center or in 
the SW region. However, this asymmetry peaks at R =2 30" and 
does not closely follow the estimated relative surface brightness 
contribution of the off-centred envelope. 

The outermost measurement from (Ventimiglia et al ■1 12010b 
of V L os = 3835 km s -1 and o- LOS = 429 km s" 1 at R = 100" 
is relevant for this discussion. These values are measured for a 
ID spectrum extracted from a region in the NGC 331 1 halo that 
is only weakly affected by the off-centred envelope; the orange- 
dashed line in Figure [T4lshows the location of this slit region on 
the image. These measurements therefore suggest that the sym- 
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R (arcsec) 

Fig. 13. Asymmetric kinematics and light distribution of the 
off-centered envelope around NGC 331 1. Lower panel: the full 
black line shows the extracted surface brightness profile from 
the VL17FORS1 V band image along P.A.= 32°, the red dashed 
line shows the Sersic profile (n = 10.5, R e = 850") of the 
MSM described in Sect. 14.21 and the green dashed line shows 
the contribution from the off-centered outer envelope along this 
P.A. The arrow indicates the dwarf galaxy HCC 027. Middle 
and upper panels: LOS velocity dispersion (cr L os, triangles) and 
velocity measureme nts (Vlos, asterisks) along P.A.= 29° from 
iRichtler et ail d2011l) (P.A.= 29° shown as R < 0). 




Fig. 14. Location of dwarf galaxies and long slit on the residual 
VLT/FORS 1 V-band image from Fig.Ql] Here the orange circles 
indicates the dwarf galaxies within 100 kpc of the NGC 3311 
center; their velocities are reported on the image. The blue circle 
shows the dwarf galaxy HCC 026, and the red diamond indicates 
HCC 007. The object to the right of HCC 007 is a foreground 
star. The yellow line centered on HCC 026 at P.A. = 142° il- 
lustrates the position and orie ntation of the long slit u sed for 
the kinematic measurements in Ventimiglia et al.l (2010) and in 
this paper. The orange-dashed section of the slit identifies the 
area where the VVos, cr measurements at R = 100" were taken 
(IVenti miglia et al. 2010); see Sect. 15.11 The red and green parts 
of the line depict the slit sections used to measure the mean ve- 
locity of the tidal stream. The FoV is 5'.1 x 2'. 8; North is up and 
East to the left. 



metric halo is hot and approximately at the galaxy systemic ve- 
locity, and that the off-centred envelope is responsible for shift- 
ing the velocities measured in the NE quadrant of the NGC 3311 
halo to more red-shifted values. We will return to this issue in 
Section[6] 

5.2. Kinematics signatures of the tidal stream around 
HCC 026 

We wish now to confirm that the stream-like substructure in the 
diffuse envelope NE of NGC 331 1 is physically associated with 
the dwarf galaxy HCC 026. The goal is to see whether we can 
detect absorption line features in the spectra at the position of 
the tidal stream, in addition to the absorption lines from the stars 
in this part of the dynamically hot halo of NGC 331 1. 

We use the deep long slit sp ectroscopic observations of the 
NGC 3311 halo carried out by IVentimiglia et al.1 (1201 Ol) with 
FORS2 on VLT-UT1. These spectra extend over the wavelength 
range 4655A - 5955A and include absorption lines from H^, 
Mgl (,L15167,5173,5184A) and Fe I (AA5270, 5328A). They 
were acquired with a 6'. 8 arcmin long-slit of width 1".6, and 
GRISM 1400+ V, giving an instrumental dispersion of 0.64 A 
pixeL 1 and a spectral resolution of cr = 90 km s~ . The 
long slit was centered on the dwarf galaxy HCC 026 at a — 
10h36m45.85s and 5 = -27d31m24.2s (J2000), with a posi- 
tion angle P. A. = 142° so that is passes through the photometric 
stream. HCC 026 itself is seen in projection onto the NGC 331 1 
halo. The slit geometry is shown in Figure[l4] Eight exposures of 



1800 sec each were taken (in total 4 hrs). In addition to the deep 
spectra, the standard G8III star HD 102070 and the spectropho- 
tometric standard star EG 274 were also observed with the same 
set-up. T he data reduction followe d standard procedures as de- 
scribed in IVentimiglia et akl (|2010). The final deep 2D spectrum 
was then used to construct radial profiles of the velocity disper- 
sion and mean LOS velocity for the halo of NGC 3311 from 
measurements at various positions along the slit. 

Here we want to additionally detect in the spectrum the ab- 
sorption lines of the stream superposed onto the main absorp- 
tion features of the hot NGC 3311 halo. We expect these sec- 
ondary absorption lines to have strengths of about 15% of those 
for the NGC 3311 halo, on the basis of the photometry carried 
out in Sect. 14.21 see Figs.[l0][TT] and[T3l If the sub-component 
projected on the outer halo is a tidal stream emerging from 
HCC 026, we expect its absorption lines to occur at significantly 
different recession velocity from those of the NGC 331 1 diffuse 
halo (Vhaio = 3921 km s ). This is based on the fact that all 
of the dwarf galaxies seen at this location, including HCC 026, 
have recession velocitie s of order 5000 km s , see Tableland 
IVentimiglia et al.1 (120111) . 

We first extract the light profile along the slit and identify 
those regions where the continuum is bright enough to provide 
suitable S/N for the kinematic measurements. We can identify 
only two such regions to the North- West of HCC 026: region I - 
7" wide (30 pixels)- and region II- 18" wide (75 pixels). These 
are indicated by the green and red slit sections in Figures H4l and 
IT5l All the flux in each of these regions is co-added to reach 
S/N - 20 in the continuum for the final one-dimensional (ID) 
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Fig. 15. Light profile along the slit centered on HCC 026, with 
the different color-shaded areas indicating the slit sections used 
for the extraction of the ID spectra of HCC 026 (blue) and the 
tidal stream (green and red). Left is SE of HCC 026, right is NW. 
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Fig. 16. ID wavelength and flux calibrated residual spectra I 
(upper panel) and II (lower panel), extracted at the position of 
the light excess. These ID spectra are obtained by subtracting 
a template model of the NGC 331 1 halo from the extracted ID 
spectra. The vertical dash-dotted lines indicate the locations of 
the main absorption features whose rest frame wavelengths are 
also reported on the figure. 



extracted spectra. We must be aware that co-adding all the signal 
from an extended portion of the slit causes a broadening of the 
absorption lines, but this does not affect our goal, which is the 
detection of a secondary absorption line component at a different 
LOS velocity from that associated with the hot halo component. 

To measure the stellar kinematics, we restrict ourselves 
to the central wavelength range 4800A < A < 5800A 



of the extracted ID sp ectrum, and use the "penalize d pixel- 
fitting" method (PPXF ICappellari & Emselleml l2004l) . In the 
PPXF method, stellar templat e stars from the MILES library 
dSanchez-Blazquez et alJl2007l) are combined to fit the ID spec- 
trum while simultaneously determining the mean LOS velocity 
and velocity dispersion from the absorption lines. From the best 
PPFX fit to the co-added spectrum we obtain values for these 
kinematic paramete rs which are cons i stent w ith the halo kine- 
matics measured by Ventimi glia et al.1 d2010l) . as well as a best- 
fit stellar template for the halo of NGC 3311. 

From the photometry carried out in Section H~2l the luminos- 
ity from the tidal stream amounts to ~ 15% of the total light in 
the region sampled by this spectrum. Attempting a direct spec- 
tral decomp osition of the stellar ha lo of NGC 3311 and the light 
excess as in lCoccato et al.l (1201 lbl) was unsuccessful. To be able 
to detect the weak kinematic signal from the excess of light, we 
therefore need to subtract the main halo contribution. We achieve 
this by taking the best-fit stellar template spectrum obtained by 
PPFX for the total extracted science spectrum, multiply it by a 
suitable fraction (0.85) and subtract it off the extracted spectrum. 
This procedure when applied to the ID spectra I and II leaves the 
residual spectra shown in Figure [TBI 

Clearly, the S/N of the residual spectra is not high enough for 
a direct pixel fitting, although the main absorption features are 
readily identified, and we must now use a different approach. We 
use the RV.FXCOR task in IRAF to identify the velocity compo- 
nents in the residual spectru m; this task implements the Fourier 
cross-correlation technique of Tonry & Davis] (1 1979b . The upper 
panels of Figure [T7] show the Fourier cross correlation functions 
computed with RV.FXCOR between the residual ID spectra I and 
II and the extracted ID spectrum of the star HD 102070. All spec- 
tra have their continuum fitted and subtracted-off, and only the 
wavelength interval 4800A < A < 5800A is used. 

In the lower panels of Fig. [17] the region centered around 
the two strongest correlation peaks is enlarged. This reveals the 
presence of two components at different velocities in both spec- 
tra: the strongest peak is found at 5054 km s , 5010 km s 
and a second peak is detected at 3931 km s _1 , 3889 km s _1 , for 
spectrum I and II respectively. The error in each velocity mea- 
surement is +55 km s . The estimated rms noise in the veloc- 
ity range 2000 - 8000 km s~'in the two lower panels of Fig. [17] 
is st 0.066 and =* 0.038. The height of both peaks in the cross 
correlation spectrum in each panel is therefore ~ 2 times the 
amplitude of the noise; because the two panels are for two in- 
dependent spectra we therefore judge these signals to be signif- 
icant. The bluer velocity peaks correspond to a component from 
the extended halo light that is still present in the residual spectra 
at an average LOS velocity of 3905 km s , see also Sect. 15.11 
The peaks at more redshifted LOS velocities indeed provide ev- 
idence for a second component at AV =* +1200 km s _1 rela- 
tive velocity with respect to the systemic velocity of NGC 331 1 
(V N GC33ii = 3800 km s _1 ), and AV =* +1100 km s _1 relative to 
the halo (V ha i = 3921 km s" 1 ). The average LOS velocity ob- 
tained by combining the measurements from spectra I and II is 
5032+38 km s -1 . F or comparison, the LOS velocity of HCC 026 
is 4946 + 4 km s~' dVentimiglia et alJl2010h . 

In summary, the stars in the stream NE of HCC 026 have 
a LOS velocity very similar to that of HCC 026 itself but 
AV =a +1200 km s~ different from the systemic velocity of 
NGC 331 1. This is a necessary condition for the stream around 
HCC 026 being a tidal stream, but could also be consistent with 
the stream stars being part of a shell which was stripped from a 
larger galaxy together with HCC 026. However , comb ining our 
kinematic result with the result of lCoccato et al.l (12.01 lal) that the 
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Fig. 17. Cross correlation functions used for kinematic measurements in two slit sections. Upper panels: Fourier cross-correlation 
functions computed between the residual (excess) spectra and the spectrum of the G dwarf star HD 102070 for spectrum I (left) and 
II (right). Lower panels: Enlarged view of the Fourier cross-correlation function peaks. The double peaks in both Fourier cross- 
correlation functions indicate the presence of two sets of absorption line features at different velocities in the residual spectra. The 
more red-shifted peaks in the Fourier cross-correlation are at velocities Vlos.i = 5054, 5010 (+55) km s , and the secondary bluer 
peaks are at V LOS ,2 = 3921, 3889 (+55) km s" 1 . 



stellar population at the NW stream position is more metal poor 
than in the surrounding outer halo, consistent with the super- 
posed stream stars having similar metallicity as the dwarf galaxy 
HCC 026 itself, makes a strong case that the stars in the stream 
were indeed tidally stripped from HCC 026. 

6. Correspondence between surface brightness 
components and kinematical structures 
measured with planetary nebulae, and the 
velocities along the HCC 007 tail 

To begin with, we summarize the results from observations of 
PN kinematics in the Hydra I cluster core, in order to see how 
they may help in unders tanding the physical str uctures identi- 
fied in the photometry. In Ventimi glia et al.l d201 ll) . we measured 
LOS velocities for 56 PNs in a 6.8' x 6'. 8 field covering the 
central (lOOkpc) 2 of the Hydra I cluster, using multi-slit imag- 
ing sp ectroscopy (MSIS, iGerhard et all 120051: 1 Arnaboldi et alj 
2007). We identified different velocity components in the PN 
LOSVD and described their spatial distributions. In brief, we 
found: 

- A broad, asymmetric central velocity component in the PNs 
LOSVD, peaked at ~ 3100 km s"'and with cr a 500 km s -1 . 
The spatial distribution of these PNs has significant substruc- 
ture, depending on the LOS velocity. It is not clear whether 
this component traces the symmetric halo around NGC 3311. 

- A minor blue-shifted velocity component, centered at ~ 
1800 kms" 1 , whose PNs have an elongated distribution 
along the North/South direction. There is no galaxy with 
similar velocity inside the 6.8' x 6'. 8 field centered on 
NGC 3311, but two such galaxies are located near the 
perimeter of the field. The association of these with the blue- 
shifted PNs is possible but unclear, however. 

- A narrow red-shifted velocity component, around Vlos — 
5000 km s -1 . These red-shifted PNs show a concentration in 



the NE quadrant of NGC 3311, and are correlated spatially 
and in LOS velocity with a group of dwarf galaxies, which 
are observed within (100 kpc) 2 of NGC 3311. 

In Figure[l8] we plot t he PNs associated with the red-shifted 
peak of the PN LOSVD dVentimiglia et al.ll201 lb on the resid- 
ual VLT/FORS1 V band image with respect to the MSM for 
NGC 331 1. This shows that a major fraction of these red-shifted 
PNs are concentrated in the NE quadrant where the offset halo 
is found. We count 10 PNs superposed on the offset envelope. 
Their average LOS velocity is V re dPNs,LOS = 5095 km s" 1 and 
CredPNs.LOS =521 km s~ . The average LOS velocity is similar 
to that of HCC 026 (V H cco26 = 4946 kms" 1 ) and to the value 
for the second component measured in Sect. 15.21 (5032 km s 
average ). With a relative offset of 1200 km s" 1 from NGC 331 1 
and a contribution of =2 15% to the surface brightness, the shift in 
the predicted mean velocity would be of order 180 km s , con- 
sistent with the measured value. However, we found in Sect. 14.31 
that the off-centered envelope contributes ~ 30% of the surface 
brightness in the NE region, and this value is consistent with the 
value indicated at the HCC 026 position in Figs. \W\ [13] Thus, 
while from their spatial location the red-shifted PNs appear to 
trace the off-centered envelope, it is not clear whether they trace 
all the stars in it. 

Therefore, from the off-centered envelope's luminosity esti- 
mated in Sect. 03] Lv,NE,env = 1.2x 10 10 (+6.0x 1O 8 )L , we can 
derive an upper limit for the luminosity specific PN number a for 
this stellar population. The measured value is a — 1/1.2X10 9 PN 
L Q l ; once we correct this number for the limitin g magnitude of 
these MSIS observations a-roT — 82 x a (see IVentimiglia 'et al.l 
d201 ll) for further details), we obtain log(axoT) = -7.18. This 
value for axoT is similar to that for old stellar pop ulations which 
are ob served in the M31 bulge and SOs galaxies (Bu zzoni et al.l 
120061) . 

For the symmetric halo of NGC 331 1 the a parameter must 
be lower. Within 25" < R < 120", we count 12 PNs from 
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Fig. 18. Positions and velocities of PNs (red diamonds, labels 
indicate the PNs' Vlos m km s~ ) associated with the red peak 
at 5000 km s" 1 in the P N LOSVD of the Hydra I cluster, from 
IVentimiglia et al.1 (1201 lh . The red peak PNs are superposed on 
the residual V-band image, which shows the substructures in the 
diffuse light in the Hydra I cluster core. The black crosses in- 
dicate the position of NGC 3311 (center) and NGC 3309 (up- 
per right), respectively. The green crosses indicate HCC 026 and 
HCC 007. The FoV is 6'.8 x 6'. 4. 



the central velocity component (IVentimiglia et al.l201 lh . If these 
trace the total luminosity of the symmetric halo, Lv.haio — 7.5 x 
10 10 (+2.0 x 10 9 ) Lq, the measured a = 1/6.2 x 10 9 PN L~ l and 
the corresponding log(o'TOT) is -7.87, a factor 5 lower than the 
value for the stellar population in the offset envelope. In fact, the 
LOSVD of the central PN component is quite axisymmetric, and 
very few PNs at the systemic velocity of NGC 3311 are found on 
the inner symmetric halo. Thus the log(a'TOT) for this component 
may be even lower. 

Superposed on the off-centered envelope, there are two tidal 
streams in the Hydra I core; they are shown in Figs. [9] and 
[TTI In Sect. 15.21 the direct measurement of Vlos of the stream 
around HCC 026 is 5032 + 38 km s , and the systemic veloc- 
ities of HCC 026 and HCC 007 are V H cc 026 = 4946 km s _1 , 
Vkcc 007 = 4830 km s respectively. Four of the PNs plotted in 
Figure[18]are consistent with being located on the NE HCC 007 
stream (one is on the polygon used to estimate the luminosity). 
These indicate a systematic velocity gradient along the stream, 
increasing from 4830 km s east of the off-centered envelope 
to a maximum of 5470 km s near the closest approach to NGC 
331 1, and decreasing again to 4890 km s near HCC 007 itself. 

Assuming that the stellar population in these streams is simi- 
lar to that in the off-set envelope, we would expect 2 + 1 PN in the 
polygon on the NE tail of HCC 007, and no PNs on the SE tail or 
on the stream around HCC 026. In Figure [18] we count one PN 
associated with the corresponding part of the HCC 007 stream, 
which is consistent with the prediction. There is not enough light 
in the stream around HCC 026 to expe ct any PN to be detected 
within the limiting magnitude of the IVentim iglia et al. 
survey. 



7. Discussion 

7. 1. The peculiar outer halo of NGC 331 1 

Outer halo or intracluster light? - Based on V-band photometry 
out to ~ 100" (~ 25 kpc), the outer halo of NGC 331 1 can be 
represented by a symmetric Sersic model with large n ^ 10, and 
an additional off-centered component with centroid shifted by 
about 50" to the NE (Sect. |4j>. The off-centered component can 
be described by a (flatter) exponential profile. However, it is not 
the extended light profile which distinguishes this central clus- 
ter galaxy from other luminous elliptical galaxies, but its very 
steeply rising velocity dispersion profile (VDP) cr(R). From a 
central value of ~ 170 km s , the VDP rises to <x =* 230 km s 
at/? = 15" - 3.7 kpc, and then on to cr = 300 - 450 km s~' a t 
R = 50" ^ 12 kpc ( IVentimiglia et al.ll2010tlRichfler et alj201 lh . 
The steep rise of the velocity dispersion profile corresponds to a 
steep increase in the enclosed mass. Within R — 20 kpc, the total 
dark matter mass inferr ed from X-ray observations is ~ 10 12 M o 
(lHavakawa et al.ll2004 . Thus NGC 331 1 is located at the center 
of the dark matter cusp of the Hydra I cluster. 

In the NW region dominated by the off-centered halo compo- 
nent, the velocity dispersion is particularly high (~ 450 km s ), 
about 60% of the galaxy velocity dispersion in the cluster core 
(IVentimiglia et alj |2010). Its mean LOS velocity is also shifted 
by ~ 100 km s" 1 with respect to NGC 331 1. Many PNs located 
on the off-centered envelope move at LOS ve locities of +1200 
km s' 1 with respect to the NGC 331 1 center dVentimiglia et al.l 
| 201ll). but also the LOS VD of the 'central PN component' of 
Ventimiglia et all (1201 ll) is quite asymmetric. The off-centered 
envelope may therefore equally well be considered as part of the 
ICL in the cluster core; because of the steeply rising VDP, distin - 
guishing outer halo from ICL is difficult (see Dolag et al. 2010). 

Dynamically hot outer halos in other BCGs? - A galaxy 
with a dynamically hot outer halo similar to NGC 3311 is 
NGC 6166, the cD galaxy in the cluster A2199 dKelson et alJ 
2002). In this system, the velocity dispersion first decreases from 
the central value of 300 km s to 200 km s within a few kilo- 
parsecs, and then steadily rises to 660 km s~ at a radius of 
60 kpc, nearly reaching the velocity dispersion of the cluster 
(o"A2i99 = 775 ± 50 km s ). This is similar but less extreme 
than the case of NGC 3311 studied here, which does not show 
the central decrease. However, such dynamically hot outer ha- 
los in BCGs are rare: e .g., for NGC 1399 in the Fornax clus- 
ter dMcNeilet alJl2010b . NGC 4874 and NGC 4989 in Coma 
dCoccato et al.ll2010al) the outer VDPs are flat, while for M87 in 
the Virg o cluster the VDP first rises out to 250" and then falls 
steeply dDohertv et al.ll2009l;lMurphv et alj|201 ll) . 

Asymmetry in X-ray observations of the Hydra I cluster 
core. - Independent evidence for asymmetries in the Hydra 
I cluster center comes from Chandra and XMM observations 
(Hayakawa et al. 2004, 2006). These authors report about an ex- 
tended emission ~ 1'.5 in the direction NE from NGC 3311 
relative to the iso-surface brightness contours further out, with 
an an gular scale of a bout 1.5', see F i g. 2 in lHavakaw a et al.l 
(120041) and Fig. 5 in lHavakawa etail (I2006I) . The X-ray sur- 
face brightness maps show a morphology reminiscent of the 
off-centered halo in the stellar light distribution. The extended 
X-ray emission NE of NGC 3311 is also the brightest among 
several high metallicity regions in the Hydra I cluster, and its as- 
sociated gas m ass is larger than that i n the compact X-ray halo 
of NGC 3311. lHavakawa et alJ (12004 interpret this high metal- 
licity region 1'.5 to the NE of NGC 3311 as the result of gas 
stripping of NGC 3311, which would imply a relative motion of 
NGC 331 1 with respect to the surrounding outer halo. It may be 
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possible that the collision with the dark matter halo of the group 
including HCC 007 could have induced such a relative motion; 
this requires further study. 

Off-centered outer envelopes in BCGs. - Separate outer 
compo nents in BCGs are frequent; the work of [Gonzalez et al.l 
(2005) for a sample of 24 clusters shows that a two component fit 
to the light profiles of BCGs provides an improved match to the 
data, and that the two p hotometric compone nts are misaligned 
in 60% of the sample. Gon zalez et al.l d2005l) reported also the 
case of Abell 1651, where the two components have different 
centers, with the outer component off-centered by about 15 kpc 
in linear distance. This is similar to the separation between the 
center of the off-centered envelope described in Sect. |4] and also 
the center of the extended X-ray emission, from the inner parts 
ofNGC3311. 

Stellar population in the outer halo of NGC 3311. - 
ICoccato et alj d201 lab measured the line strength indices for 
the stellar population in the NGC 3311 halo along a 6'. 8 slit 
centred on HCC 026 and with P.A.=142°. These measurements 
show the effect of the stars in the HCC 026 tail: in this region, 
the metallicity is lower, [Fe/H]= -0.73 + 0.06, compared to 
the average value for the outer halo away from the tail, where 
[Fe/H]haio - 0.34 + 0.05. The metallicity in the tail region is 
consistent with a ~ 20% contribution of stars with the metal- 
licity of HCC 026. Unfortunately, in the halo region away from 
the HCC 026 tail, these data do not have sufficient S/N to test 
whether the stars in the offset envelope are a different stellar 
population from those in the symmetric halo. 

7.2. Tidal tails and streams around NGC 331 1 

Comparison of tails in Hydra I with those around other galax- 
ies. - It is relevant to compare the luminosities of the streams 
in the Hydra I core with those o bserved in the halos of bright 
ellipticals in the Virgo cluster by iJanowiecki et all (1201 Ol) . The 
HCC 026 stream identified here with total luminosity slightly 
below 10 9 L o ,v is significantly brighter than any of the Virgo 
streams, and the HCC 007 stream is several times more lumi- 
nous still. Also, the tidal streams seen around the Virgo cluster 
ellipticals have much fainter peak surface brightness (brightest: 
pv = 25.7 mag arcsec~ 2 ; average: py = 27.6 mag arcsec~ 2 ) 
than those measured in the Hydra I streams, which are py = 
24.8 mag arcsec~ 2 for the HCC 026 stream and py = 24.4 
mag arcsec~ 2 for the HCC 007 stream. 

Streams with as low surface brightness as those observed in 
Virgo would not be detectable in our Hydra images, thus could 
easily be abundant in this region. In fact, the two small objects 
approximately near the middle of the polygon on the NE stream 
of HCC 007 in Fig. Qj] may be connected by a small stream 
within the HCC 007 tail, as suggested especially by the right 
panel of Fig. [9] This small stream, if confirmed, would have a 
luminosity more similar to the Virgo streams, albeit a higher sur- 
face brightness. 

We can expand this comparison to the streams a round 
bright local elliptical galaxies analyzed bv lTal et all (|2009). The 
Hydra I streams are detected at radii larger than ~ 23 kpc 
from the center of NGC 3311, and the HCC 007 tail reaches 
at least 110 kpc in total projected length, see the measurement 
in Sect. 14. II These numbers for the Hydra I streams are consis- 
tent w ith the av e rage r adius of occurrence of the streams stud- 
ied by iTal et al. | (120091). and t he linear extent of some of them. 
Unfortunately,|TaleLalJ d2009) do not give luminosities for their 
streams. 



The streams are tidal tails. - The streams emerging from 
HCC 026 and HCC 007 can be interpreted as the result of tidal 
disruption. In the case of HCC 007, the sheer size of the stream 
(the two tails together span ~ 110 kpc) is difficult to explain 
by any other model. For HCC 026, the similar LOS velocity of 
the stream with the galaxy HCC 026 could still allow a model in 
which the stream stars are part of a shell which was stripped from 
a larger galaxy together with HCC 026. However. ICoccato et al.l 
(201 la) showed that the metallicity obtained with a single stellar 
population model in the stream region is lower than in the dif- 
fuse stellar halo of NGC 3311 away from the stream, and that 
this lower metallicity value is consistent with a composite pop- 
ulation of stars from HCC 026 and stars from the stellar halo of 
NGC 3311 in approximately the proportion implied by the sur- 
face brightnesses of both components. This makes a strong case 
that the stars in the HCC 026 stream were indeed tidally stripped 
from HCC 026. It also shows that a fruitful way to investigate the 
origin of the diffuse light features in Hydra I and elsewhere is to 
measure stellar population parameters from deep spectroscopy. 

Given their redshifted velocities, both the HCC 007 and 
HCC 026 streams would likely be located now on the distant 
side of the Hydra I cluster core, with their stars being on slightly 
different orbits than the galaxies they were once bound to. 

Comparison with simulated tails. - The overall morphology 
of the two streams in the Hydra I cluster core is si milar to those 
of tid al streams in the simulation investigated by Rudick et al. 
(2009) to study the formation of intracluster light. In particular, 
the tidal streams in their case Gl reproduce one characteristic 
property of the streams discovered here, especially for HCC 007, 
that one side of the st ream is significantl y brighter than the other. 
In the simulation of iRudick et alj d2009l) . the ~ 200 kpc long 
tidal arms in Gl formed as a consequence of the interaction of 
a disk galaxy with the cluster cD, with a pericenter distance of 
~ 100 kpc. Most of the stars in the tidal streams in this system 
were unbound from the disk galaxy shortly after the pericentre 
passage, and once formed, the streams decay in ~ 1.5 orbital 
times. The fact that tidal streams are visible around HCC 026 and 
HCC 007 thus indicates that both galaxies have recently passed 
the pericenters of their orbits. 

7.3. A group of galaxies in disruption in the Hydra I cluster 
core 

All the galaxies in d isruption. - As already reported by 
IVenti miglia et al] (1201 ll) . there are no galaxies in the Hydra I 
core at cluster-centric radii (< 100 kpc) with velocities around 
the cluster systemic velocity. However, several dwarf galaxies at 
high velocities V LO s > 4500 km s'and the SO galaxy HCC 007 
are located in this region. In Table [3] we list the sky coordinates, 
apparent total V-band magni t udes, and LOS velocities for these 
galaxies from Misg eld et al.l d2008l) . These galaxies are part of 
a well-defined cluster substructure bot h in velocity and spatial 
distribution, as already commented by Ventimig lia et al.l (1201 ll) 
and shown in Figure[l4] For two of these galaxies, HCC 026 and 
HCC 007 we have been able to detect tidal tails, showing that 
they are being disrupted by the tidal field near the cluster center. 
For the others, the tidal effects may not yet be strong enough to 
have lead to detectable tidal tails, or their tails have dispersed 
below the detection limit. Independent of this, the observed tidal 
streams are the consequence of the recent infall of a cluster sub- 
structure containing the galaxies at Vlos > 4500 km s . 

Where have the central galaxies gone? - We now return 
to the observed absence of galaxies at cluster-systemic veloc- 
ity in the central 100 kpc around NGC 3311 (a similar result 
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Table 3. Dwarf galaxies in the NGC 3311 field. Magnitudes and line-of-sight velocities are from Misge ldet al.1 (12008). (*) value 
from Ventimiglia (201 1), Ph.D. Thesis. 



was found in the NGC 5044 group, iMendel et ail (120091) 1. A 
plausible explanation is that such galaxies are no longer seen 
in the central region of the cluster because they were all dis- 
rupted in the past during close encounters with the central galaxy 
NGC 3311 and with the dark m atter cusp at the cluster center 
dFaltenbacher & Mathewsll2005l) . In this case, their former stars 
would now be part of the diffuse stellar component in the Hydra I 
core, including the halo around NGC 3311. It is interesting to 
ask whether most of the stars in the outer halo could have origi- 
nated from disruption of small galaxies in the way that seems to 
be currently on-going, or whether the majority of the halo stars 
come through a different channel, the tidal disruption of the ha- 
los of massive elliptical galaxies prior to merging with the clus- 
ter cD dMurante etaDl2007t IPuchwein et alll2010l) . Again, anal- 
ysis of the stellar population properties of dwarf galaxies and 
of the outer halo stars is imp ortant for answering this question 
dCoccato et alJl2010bll20TTal) . 

Phase-mixing and formation of ICL from the tidally dis- 
solved stars. - Assuming that the interaction with NGC 3311 
and the Hydra I cluster core ultimately disrupts the dwarf 
galaxies seen close to NGC 3311 completely, all of their stars 
would phase-mix and distribute in the central cluster poten- 
tial. Converting the uni-directional velocity of their orbits of 
now ~ 1200 km s^relative to NGC 3311 to 3o-;L , the result- 

ing cr eq - 700 km s . This would be somewhat increased by 
adding the contribution from the unknown transverse velocities. 
On the other hand, some of the kinetic energy of these stars 
would be converted to potential energy during the phase-mixing, 
lowering the final cr eq . For comparison, the velocity dispersion 
of the NGC 331 1 halo near the current position of HCC 026 is 
Chaio — 400 km s . This suggests that the stars from the galax- 
ies that are presently tidally disrupted will end up at larger radii 
than their current position, and with somewhat larger dispersion 
than the observed cr ha i a 400 km s . I.e., they will end up in 
the outermost halo and intracluster light around NGC 331 1 and 
in the Hydra I cluster core. 



8. Summary and conclusions 

In this work we extend our investigation of the properties and 
origin of the diffuse light in the Hydra I cluster. Combining 
surface photometry with kinematic information from long slit 
and planetary nebula (PN) data, we find an off-centered, diffuse 
outer halo around the central galaxy NGC 3311, and show that 
at least two galaxies are currently disrupted in the cluster core, 
adding their stars to the outer halo and intracluster light around 
NGC 3311. More specifically, our results are as follows: 

Structural parameters are derived for the two giant elliptical 
galaxies in the cluster core, NGC 3309 and NGC 331 1, using V- 



band imaging data obtained at the ESO/MPI 2.2m telescope and 
archival VLT/FORS1 V band data. While the light distribution 
of NGC 3309 is reproduced by a single Sersic profile, that of 
NGC 331 1 is characterized by several components. Outside the 
nuclear regions which are affected by a dust lane and by bright 
luminous knots, the bright regions within 30" follow an R 1 ^ 4 law. 
In the deep V band data, these regions together with the symmet- 
ric part of the outer halo can be described by a Sersic law with 
n = 10. 

The residual image, obtained after subtracting the two- 
dimensional model of the two bright galaxies, shows an addi- 
tional extended envelope centered at ~ 80" to the North-East 
of NGC 3311. This off-centered envelope is approximately de- 
scribed by a exponential profile and contains Lv,NE,env = 1.2 X 
10 10 (+6.0 x 1O 8 )L . This corresponds to ~ 50% of the lumi- 
nosity of the symmetric halo in the same region, and ~ 15% of 
the luminosity of the entire symmetric halo in the radial range 
25" - 120". 

Furthermore, the diffuse light in the Hydra I core harbors 
two tidal streams emerging from the dwarf galaxy HCC 026 and 
the SO galaxy HCC 007 (see Fig. [TTJ. The total luminosity in 
the NW part of the HCC 026 stream is L v ,nw,hcc 026 = 4.8 x 
10 8 (+8xl0 7 ) Lq, with an average surface brightness jjy = 24.8 + 
0.2 mag arcsec 2 , ~ 15% of that in the extended envelope and 
symmetric halo at the stream position. 

Analysis of a deep spectrum in this region indicates that 
the NW stream has a LOS-velocity similar to HCC 026 itself 
(~ 5000 km s , see Sect. |5}, and that it consist s of stars with 
similar metallicities as the stars of HCC 026 (Coccato et al. 
1201 lab . These results favour the interpretation that the HCC 026 
stream has been tidally dissolved from HCC 026, over one where 
both the dwarf galaxy and the stream were dissolved from a 
larger galaxy. The luminosity of the combined NW and SE parts 
of the stream is several times the current luminosity of HCC 026, 
^v,hcc 026 = 1.5 X 1O 8 L , so that this dwarf galaxy has by now 
been mostly dissolved by the tidal field. Both HCC 026 and its 
tidal tails have a relative LOS velocity of ~ 1200 km s with 
respect to NGC 331 1 (V N GC33ii - 3800 km s" 1 ). 

The second tidal stream around the SO galaxy HCC 007 
extends over at least ~ 110 kpc on both sides of this galaxy. 
The tail is fairly thick and is brighter on the NE side where 
fiv = 24.4 + 0.5 mag arcsec~ 2 . On the NW side it is about one 
magnitude fainter and less certain. We measured a luminosity of 
iv,NE,HCC 007 = 2.9 x 10 9 (+5 x 10 8 ) L in the brightest region of 
the NE tail. The NE tail appears to join the outer part of the off- 
centered halo east of NGC 33 1 1 ; if so, its total luminosity could 
be significantly higher. For comparison, the total current lumi- 
nosity of HCC 007 is ~ 4.7 x 1O 9 L . This galaxy has therefore 
lost of order 50% of its stars. 
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Four PNs are found on the NE HCC 007 tail which allow 
a preliminary measurement of its LOS velocity. The largest PN 
velocity (5470 km s _1 ) is found SE of NGC 3311 in a region 
dominated by tail stars. One PN is close to HCC 007 itself and 
has a velocity within 60 km s~ : of the galaxy's systemic velocity. 

The morphologies of these streams are similar to those of 
tidal streams generated when galaxies are disrupted on highly 
radial orbits through the cluster cent er, shortly after the time of 
closest approach with the cluster cD dRudick et al.ll2009l) . 

Superposed on the off-centered outer halo of NGC 33 1 1, we 
also find a number of PNs from the same sub-component of the 
PN LOSVD centered at 5000 km s , as well as a group of dwarf 
galaxies with similar velocities. This suggests a physical associ- 
ation between the off-centered halo, the two tidal streams, and an 
entire group of galaxies at about 5000 km s systemic velocity, 
including HCC 026, several other dwarf galaxies, and HCC 007, 
which are currently falling through the core of the Hydra I clus- 
ter and have already been partially disrupted. 

What fraction of the off-centered halo is due to stars dis- 
solved from these galaxies is, however, not constrained well 
by our measurements. The deep spectra show only a modest 
shift in mean LOS velocity in the NE halo, and the photomet- 
ric, kinematic and stellar population analysis of the HCC 026 
tail suggests that the component at 5000 km s LOS velocity 
contributes only a fraction of the inferred off-centered halo at 
this location. Furthermore, a number of PNs with moderately 
to highly blue-shifted ve locities are also seen su perposed on the 
outer halo of NGC 3311 dVentimiglia et al.l201 ll) . To resolve this 
requires further kinematics and stellar population analysis with 
deep spectroscopy at different positions in the halo. 

This work provides a vivid example of how strong tidal 
forces cause morphological transformations of galaxies falling 
through the cores of galaxy clusters. It also shows that stars are 
currently being added to the diffuse light in the Hydra I cluster 
core. The stars in the streams around HCC 026 and HCC 007 
must have been unbound from their parent galaxies during the 
on-going close passage through the high-density center of the 
cluster, and are now on slightly different orbits from their once 
parent galaxies. Ultimately they will phase-mix in the potential 
and add to the diffuse outer envelope of NGC 3311 and to the 
intracluster light in the cluster core. 

The on-going accretion of a group of galaxies will add about 
several 10 9 L o .y to the outer halo and intracluster light around 
NGC 3311, about 5% of the total light in the NGC 3311 symmet- 
ric halo. Tidal dissolution of small galaxies is thus an important 
channel for creating diffuse light in clusters. It is unlikely, how- 
ever, that this is the main mechanism for the orig in of the sym- 
metric halo around NGC 33 11. Both simulations (Murante et al.l 
l2007t IPuchwei n et al. 2010) and studies of diffuse light in clus- 
ters like Coma an d the star pile cluster dGerhard et al.l l2007l 
ISalinas et al.ll201 lb indicate that the disruption and merger of a 
giant galaxy is the more efficient route to creating large amounts 
of diffuse light. Deep spectroscopy to measure Lick indices and 
stellar abundances in different regions of the offset envelope and 
halo of NGC 3311 will be important to address this issue further. 
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